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j.^ . Abstract 

O 

(T^ ■ We describe a simple way of incorporating fluctuations of the Hubble scale during the horizon exit 
\f^ . of scalar perturbations into the 5N formalism. The dominant effect comes from the dependence of 
O ' the Hubble scale on low-frequency modes of the infiaton. This modifies the coefficient of the log- 
enhanced term appearing in the curvature spectrum at second order in field fiuctuations. With this 
modification, the relevant coefficient turns out to be proportional to the second derivative of the 
tree-level spectrum with respect to the infiaton ip at horizon exit. A logarithm with precisely the 
same coefficient appears in a calculation of the log-enhancement of the curvature spectrum based 
5^ \ purely on the geometry of the reheating surface. We take this agreement as strong support for the 
proposed implementation of the 5N formalism. Moreover, our analysis makes it apparent that the 
log-enhancement of the infiationary power-spectrum is indeed physical if this quantity is defined 
using a global coordinate system on the reheating surface (or any other post-infiationary surface 
of constant energy density). However, it can be avoided by defining the spectrum using invariant 
distances on this surface. 



O 
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1 Introduction 

It is well known that curvature perturbations created during cosmological inflation [Tl|2] are affected 
by infrared (IR) divergences [3|l5j (see [6j for some recent discussions). These divergences are closely 
related to the familiar divergence of the scalar- field correlator in de Sitter space [7]. While, at leading 
order, the divergence can be absorbed into the definition of the background and is hence unobservable, 
higher orders in the curvature perturbation lead to corresponding log-enhanced corrections to the 
power spectrum. Since the IR cutoff appearing in these logarithms is provided by the size of the 
observed universe (rather than by, for example, the size of the universe created by inflation), these 
logarithms are, however, not particularly large in practice j5]. Nevertheless, it is conceivable that 
the power spectrum is measured by a very late observer, who has access to the entire region of the 
universe created in 'our' inflationary patchl^ For such a 'late' observer, the infrared logarithms can 
be extremely large and it is an interesting question of principle how to achieve consistency between 
his and our (i.e. the 'early' observer's) measurement of the power spectrum. 

We approach these issues starting from the SN formalism [S1[S] o It turns out that the problem 
can be resolved rather easily if a simple modification of the 6N formalism, which takes fluctuations of 
the Hubble scale during slow- roll inflation into account, is implemented. In essence, this modification 
consists of treating the Hubble scale, which defines the normalization of the scalar-field correlator, 
as a function of the perturbed background value of ip relevant at the time of horizon exit of a given 
mode. Since this background perturbation depends only on modes with a smaller wave number, our 
proposal can be implemented in an unambiguous and straightforward way. 

The Hubble-scale fluctuations discussed above lead to slow-roll suppressed but log-enhanced 
contributions to the power-spectrum, similar to the familiar 'c-loop' effects of the SN formalism. 
Combining both these contributions, the coefficient of the first log-enhanced correction to the power 
spectrum takes a very simple form. It is essentially given by the second derivative in (f of the 
leading order power spectrum, N^{(f)'^H{(py/(27Ty. Here the argument (p is the value of the infiaton 
corresponding to the horizon exit of the wave number k under consideration. 

The physical significance of this proposal can be understood as follows: Consider the reheating 
surface (or any other surface of constant energy density after the end of infiation). Disregarding vector 
and tensor modes, the geometry of this surface, parameterized by coordinates y = {y^ , y'^ , y^) , can 
be characterized by a single scalar function C(^) known as the curvature perturbation. The power 
spectrum may be defined as the logarithmic derivative of the correlation function {C{x)({x + y)) with 
respect to y = \y\. Alternatively, a closely related spectrum based on a fixed invariant distance s 
between pairs of points in the correlator can be defined. This latter power spectrum is an entirely 
local quantity and its expectation value does not depend on the size of the region over which the 
correlator is measured (i.e. on the 'age' of the observer). The two spectra, which we denote by V(^{y) 
and Vi^is), are related by 'P^(y) = (P^(ye'')), where C is a coarse-grained value of ( relevant for the 
distance-measurement between a given pair of points. In the relation between the two spectra, the 
averaging is over the potentially large observed region, with large expectation values of ( and C^ 
in the case of a late observer. This leads to a log-enhancement of V(^{y) which is easily seen to be 



^ Of course, this possibility is in practice limited by the presently observed dark energy or cosmological constant. 
^ The physical importance of higher-order terms in the 5N formalism was first appreciated in [S] . 



precisely the log-enhancement found earlier in (our implementation of) the SN formalism. 

On the one hand, this agreement provides support for the implementation of the 6N formalism 
that we advocate. On the other hand, it makes the 'physical reality' of large logs from IR divergences 
particularly clear: The log-enhancement arises due to the use of global coordinates in a very large 
region, where these deviate significantly from the invariant distance. It can be avoided if one measures 
the power spectrum 'Pf(s), which is defined using a two-point correlation function based on the 
invariant distance between each pair of points appearing in the spatial average. 



While our work was being finalized, Ref. [10] appeared, which overlaps with part of our analysis. 
We will comment on this in more detail at the end of Sec. El 



2 Hubble scale fluctuations in the 6N formalism 

During inflation, the amplitude of scalar field fluctuations is controlled by the Hubble parameter H. 
The latter is usually evaluated at the value ip(tk) of the classical homogenous solution ip, where tk 
is the time of horizon exit of the mode k: H = H{ip{tk)). However, this approach does not account 
for the fact that perturbations with wavelength larger than k~^ have already left the horizon. These 
perturbations modify the value of the scalar field relevant for the mode k and need to be taken into 
account. 

We do so by writing the scalar field perturbation in a quasi-de-Sitter background as 

7^77= H {^{tu) + 5^{x)) a^, (1) 

where 

^^^^) = (^^77^3 ^^^^'')) ^r, (2) 



(27r)3 72P 

where ag is a normalized Gaussian random variable q In writing eq. ([T]) and in the rest of the paper, 
we only include the leading order contributions to the fluctuations and neglect slow-roll corrections 
that are not enhanced by potentially large logarithms. The 6ip{x) contribution in the argument of the 
Hubble parameter accounts for the backreaction of long wavelength modes on the scalar fluctuations 
6ip{x) generated inside the horizon. 



In order to analyze their effect, we expand the Hubble parameter up to second order in 6ip, 
where we have introduced the Gaussian field 



M.)-/(0;e-\, 



«-^"- « 



^It satisfies the relation {a-^a-p) = {2tt)'^ 5'^'^\k + p) . 



We use the abbreviations Ai = H^p/H and A2 = H^^/{2H) for coefficients involving the ffist and 
second derivative of H with respect to ip. Note that as a result of the corrections proportional to 
Ai and A2, the scalar fluctuation in Fourier space, S(fj^, are not Gaussian at Hubble exit. A similar 
idea was considered to second order in the scalar field perturbations, in the context of the tree-level 
bispectrum in [11]. Note that we neglect the quantum mechanically generated non-Gaussianity of 
the fields which is present at horizon exit [T2] . 

We now turn to the curvature fluctuation ( evaluated on a surface of uniform energy density (for 
our purposes, the reheating surface). The 6N formalism calculates ( starting from scalar fluctuations 
on a flat surface, 



C(x) = N{^ + 6^{x)) - N{^) = N^6^{x) + - N^^6^{xY + -N^^^6^{xf +■■■ . (5) 



-N^^6ip{xy + -. 

Here N^^, N^p^ are derivatives of the number of e-folds A^ as a function of (f, evaluated on the classical 
background trajectory. The resulting two point function for curvature perturbations in Fourier space 
is given by 

iCQ) = K(5^,5^,) + N,N,,{5^,{5^')^) + ^{ {5^\ (V),) + ^^^(%.(V),). (6) 

Note that correlators of odd powers of 5ip are in general non-zero since (Jy?^ is not Gaussian. However, 
one may check that nevertheless {5^^) = 0. We have not set < C, >= 0, though making this choice 
would not change in an essential way our results (we will discuss this in more detail later). 

We now define the spectrum V^ through (x^Xp) = i2'K)^6^{k + 'p)2TT'^V^{k)/k^. This implies 
that Vy^ik) = (if/27r)^. The spectrum of curvature perturbations P^ is defined in an analogous way. 
It follows from a straightforward evaluation of eq. (E]) (see Appendix A for details) that 

V^k) = V^ [Nl + V^ HkL) [{Aj + 2A2) Nl + AA^N^N^^ + Nl^ + N^N^^^] ] 

Here it is crucial that the IR cutoff L is set by the size of the region in which the correlator is 
measured. As mentioned before, we take it for granted that perturbations on even larger scales can 
be absorbed in the background and are hence unobservable. In other words, the classical evolution 
starts with the horizon exit of modes of scale k ~ 1/L. Our main interest is the dependence of the 
spectrum on the size L of the observed region. 

We can express the second derivative along the scalar field, appearing in equation ([7]), in terms 
of slow-roll parameters. We define as usual e = V^/ {^V"^), r] = V^p^/V, ^'^ = V^V^^.^/V'^, and 
we set Mp = Svr. Using the well known relations n^ — 1 = 27] — 6e for the spectral index, and 
a^ = 16e?7 — 24e^ —2^^ for its running, i.e. a^ = dn^/rflnA; (all quantities evaluated at Hubble exit) 
we find 

^ (A^X) = ^[(V- 2e) (nc - 1) + (n, - 1)^ + a,] (8) 

This quantity vanishes in the case of scale invariance, i.e. when n^ = 1 and a^ = 0. In this case 
logarithmic corrections to the power spectrum drop out, at leading order in the slow-roll expansion. 
We will discuss an interpretation of this result in the next section. 



At our level of accuracy, we can rewrite eq. ([7]) as 



-Pdk) = V^Nl + {^%, i^ (iVX) , (9) 



where 

H/k 



(V)iA= / (^^=^xMA:L) (10) 



is the expectation value of Sip^ measured on a length scale 1/k. We see that the log-enhanced cor- 
rection takes the suggestive form of the second term in a Taylor expansion. The physical meaning of 
this structure will be clarified below. 

We note that the integration in eq. fITU]) was performed assuming that (p has an exactly scale 
invariant spectrum. More generally, if the power spectrum of (p has a constant spectral index n^p, the 
\n.{kL) term above should be replaced by 



ln(A;L) ^ —^ (l - (/tL)-^"^"'^) . (11) 



However, provided that we do not consider an exponentially large range of scales and the field 
fluctuations are reasonably close to scale invariant, then \n^ — 1| ln(/cL) <^ 1, and the more general 
result above is well approximated by \n{kL). Since the running of the spectral index is expected to 
be of order (n^ — 1)^, further corrections associated with this running are suppressed by a similar 
argument. 



3 Infrared divergences from the geometry of the reheating 
surface 

We now provide a physical interpretation for the log-enhanced correction to the power spectrum 
given in eq. ([7]). An observer, in order to make a measurement, specifies a coordinate system on the 
reheating surface (or on any other surface of constant energy density after the end of inflation). This 
can be done as follows. Choose a coordinate system where slices of constant time t have uniform 
energy density. Neglecting vector and tensor modes, we write the metric as 

dsl = -dt^ + a^{t) e^'^^^) 5ij dx'dx^ . (12) 

We use conventions where reheating occurs at time t/, at which the homogeneous scale factor a{tf) = 
1. This leads to the following metric on the three-dimensional reheating surface: 

dsj = e^^^^) 6ij dx'dx^ . (13) 

Important consequences for the power spectrum of ( can be derived just from the geometry of 
the reheating surface specified above. To see this, let us define the power spectrum as the logarithmic 



derivative of the two point function, 

vdy) = -7^\{{m-as + y)f) = --^{mas+y)). (14) 

amy z amy 

where y is the length of a coordinate vector on the reheating surface, y'^ = {y^Y + iv'^Y + in'^Y- It is 
not difficult to see that eq. f lT^ gives the correct value for the power spectrum associated with the 
fluctuations of a slowly-rolling scalar in quasi-de Sitter background (see Appendix B). 

Alternatively, we can define the power spectrum averaging over pairs of points separated by a 
fixed invariant distance s, 

Vds) = ^^ {{m - as + y{s)f) . (15) 

Here y{s) denotes a coordinate vector with invariant length s. We have introduced a tilde to dis- 
tinguish this spectrum from eq. f lT4|) . It is clear that long- wavelength background fiuctuations of (, 
which affect only the parameterization but not the physics of any local patch of the reheating surface, 
are irrelevant for this new spectrum. In other words, 'Pf(s) is a purely local quantity which, by its 
very definition, cannot depend on the size L of the observed region. 

Defining ( as the average value of ( characteristic of a small region containing a particular pair 
of points X and x + y{s), we can write 

^<w-d^K(^*''>"^*''+ '''■'))■ (!«' 

Here e is a coordinate unit vector. We see that the two spectra differ only by the selection of pairs 
of points in the averaging procedure. In one case, pairs with fixed coordinate distance y, in the other 
case pairs with fixed invariant distance 

s = eS (17) 

are chosen. Hence we can express Vc in terms of the locally defined spectrum Vr through 



Vdy) = -^ iCix) C{x + y)) = -^ (C(f) C {x + e{ye^)e-^)) = (Vdye^)) . (18) 



Here the second equality holds because logarithmic derivatives in y and s agree. The third equality 
relies on the fact that "P^ is a local quantity, as explained before. The averaging in the last term 
remains non-trivial because of the variation of C over the large patch of size L on which P^ was 
originally defined. 

Expanding to second order in (, we find 

vdv) = VM + iO^/dy) + i (?) ^.^dy) ■ (w) 

This quantifies the deviation of the observer-patch-dependent spectrum P^ from the locally defined 
spectrum P^ for large patches. It also clarifies the origin of the large logarithms of L found in the 6N 



formalism. They arise because the expectation values of ( and C^ grow logarithmically with L. We 
note that, if the power spectrum P^ is scale independent, the two last terms in eq. (fT9|l vanish. This is 
compatible with eqs. ([7]), ([H]) and has a simple intuitive reason: The log-enhancement occurs because 
in regions with a large background value of ( the scale of a given mode is effectively misidentified. 
However, in the scale-invariant case, such a misidentification has no consequences. 

We now compare eqs. (IT9l) and ([7]) quantitatively. To achieve this, we can rewrite derivatives in 
\ny in terms of derivatives in (f using dln{y)/d(f = N^ (alternatively we can directly use the results 
of Appendix B) and express expectation values of ( and ("^ in terms of (y9- correlators using eq. ([5]). 
Equivalently, we can observe that eq. f ll9p is simply a second order Taylor expansion around the 
classical trajectory. Since there is an unambiguous functional relation between fluctuations in ip and 
in C (the latter being equivalent to In y) , we immediately conclude that 






Mv) = 'Pdy)+o-!-f i'>v'} = My) + o ^x HlM -t4 ■ (20) 



Here (^v^^) is the zero-momentum 6ip correlator with UV cutoff y and IR cutoff L, in analogy to 
(C^) above. Furthermore, we used the fact that the expectation value of Sip vanishes. Thus, we 
have achieved complete agreement with the IR divergences or, more precisely, with the log-enhanced 
corrections which arise in the 6N formalism. 

To summarize, we have obtained a simple physical interpretation for the logarithmic contributions 
to the power spectrum that were determined in the previous section: The log-enhancement arises 
due to the use of global coordinates in a very large region, where these deviate significantly from 
the invariant distance. It can be completely avoided if the power spectrum "P^ is considered, which is 
based on the two-point correlation function defined in terms of the invariant distance between each 
pair of points appearing in the spatial average. 

As mentioned in the Introduction, a discussion closely related to the present section appeared 
in [To] while this paper was being finalized. In particular, using an argument slightly different from 
ours, a log-enhancement arising from the large-scale geometry of the reheating surface was derived. 
Our findings can be brought into agreement with [10] if we drop the term ~ {() in eq. (1191) and set 
a^ = dn(^/d\n k to zero. This corresponds to keeping just the (n^ — 1)^ term in our eq. (jHD- Indeed, the 
term proportional to {() disappears if we rescale our coordinates on the reheating surface according 
to y — )■ ye~''. This can be directly checked from eq. flT9|) and corresponds to the fact that, in these 
new coordinates, the average of ( vanishes. While this is certainly a rather natural coordinate choice, 
it does not serve our purpose of comparing with the 6N formalism: Indeed, in this new coordinate 
system the value C = does not any more correspond to the endpoint of the classical trajectory. 
Nevertheless, setting questions related to the 6N formalism aside, the term ~ {() appears just to 
be matter of different conventions between [10] and our discussion. By contrast, we believe that our 
term ~ a^ is real and, in general, no more slow-roll suppressed than the (n^ — 1)^ term. 

The agreement of the findings of this section (and of [lOj, subject to the caveats mentioned above) 
with the previously discussed calculation in the 6N formalism is a non-trivial result. In particular, [TU] 
argue that any attempt to understand the IR divergences using the 6N approach is incomplete due 
to certain divergences of the Sip correlation function which are already present at horizon crossing. 
One way of interpreting our modification of the 6N formalism is by saying that, via Eq. ([3]), we 



include such divergences in the formahsmo It thus appears that our purely classical modification to 
the 6N formalism is sufficient to fully capture the leading logarithmic behaviour. 

4 Discussion 

We have studied IR divergences during inflation using both the 6N formalism and a simple, phe- 
nomenological approach based just on the geometry of the reheating surface. By implementing a 
simple modiflcation of the SN formalism, we took into account the effect of modes that left the 
horizon long before the scales we are observing on the Hubble scale. Including this effect provides 
new log-enhanced contributions to the power spectrum, at the same order in H and slow-roll as 
the standard classical loop corrections. We found that the combination of all contributions can be 
assembled in an elegant formula, in which the log-enhanced contributions are weighted by the second 
derivative of the tree level power spectrum, with respect to the inffaton field. 

This result can be understood intuitively by considering two power spectra: One is defined locally 
on the surface of reheating, using invariant distances to define the correlator. The other is based on 
the coordinate distance on this surface and depends on global features of this surface, in particular on 
long-wavelength modes. When expressed in terms of the local spectrum, this latter, global spectrum 
exhibits an IR divergence associated to the size of the region on which it is measured. It is, in 
fact, this latter spectrum that is calculated in the 6N formalism and the log-divergence found in 
both approaches is precisely the same. This provides strong support for the modification of the 6N 
formalism we propose. In the case of an exactly scale invariant spectrum, the IR logarithms are 
absent. For an observer dealing with a scale-dependent spectrum and having a very large region 
available for his measurement, the use of the local spectrum, which is not affected by our IR effects, 
appears to be clearly favored. 

In single-field, slow-roll infiation the coefficient multiplying the logarithm is heavily suppressed. 
While it is usually also suppressed in multi-field models, this statement does not hold any more in 
complete generality. In fact, in some special cases the coefficient could be large enough to compensate 
for the power spectrum suppression found in single-field models [TStfl^. It would therefore clearly be 
of interest to extend our results to multi-field infiation. In these models, which contain isocurvature 
perturbations during infiation, the Hubble scale will depend on several fields, even if the primordial 
curvature perturbation depends on only one field. In cases where the loop correction is not small, 
it has been argued they could give an observable contribution to ( through non-Gaussianity [9|[l5]. 
in particular through a special type of scale dependence of the bispectrum non-linearity parameter 
[T3l[ll]. For a discussion of the scale dependence of the tree level bispectrum see [16]. It would thus 
be interesting to study the effect of our loop corrections on the bispectrum (and higher order ra-point 
correlators) . 



^ We would like to thank Martin Sloth for a conversation helping us to correct our comments on this issue, in the 
first version of the paper. 
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Appendix A: Correlator of the curvature perturbation 

In this appendix, we present the calculation relating the correlator of curvature perturbations ( to 
correlators of the Gaussian random field x- The terms which need to be calculated are 

In the previous equations, expressions like (^v^^)^ indicate convolutions. The terms in the last two 
lines contains correlators between convolutions, that provide (see the first reference in [15] for more 
details) 

( (V),. (V),^) = 4P,ln(fcL) ,{xjiXp) 

Using the expansion of eq. ([3]), written in momentum space, the 2-point correlator of the scalar field 
Sip is given by 



[^V^k^Vp) = iXkXp) 



,^+(^^+'^^^t«.(^^^^'^^-^ 



(23) 



{XkXp)[l + {Al + 2A,)V^HkL)] 



Note that, after applying Wick's theorem, only one term survives the conditions l,m <^ k. The 
remaining correlator ( Sipj: {6ip'^)^), appearing in the second term of eq. fl2T]) . is of order three in 6ip. 
Therefore, only terms at next-to-leading order, in the expansion of eq. (|3]), contribute to it. One finds 

f d^q f dH 

. f d^Q f dH , , ^ ' 

+2^1 / ^^^^i / 77;z^^XkXiXa--iXp-q)- 



, (27r)3 7,<<^ (27r)3 

Decomposing the four point functions appearing in the integrands by means of Wick's theorem, and 
taking care of the conditions on the size of /, one finds 

( 5^,^ ( V),^) = ^A{P^ HkL) ( XkXp) ■ (25) 

9 



Putting all terms together, the final result for the spectrum of curvature perturbations reads 

V^{k) = V^ {iVj + V^ HkL) [{Aj + 2A2) Nl + AA^N^N^^ + iV^^ + N^N^^^] } 

1^ /,.,^ X (26) 

2 V 



r,Nl + v,HkL)-—{Nlv,). 



Appendix B: Power spectrum in coordinate space 

In this appendix, we discuss properties of the power spectrum in coordinate space, defined in eq. ([Hj) 
as 

n(l/) = -^(C(5^)C(x + y)) • (27) 

The power spectrum in coordinate space is no more difficult to handle than the usual power spectrum 
in momentum space. Moreover, as discussed in the main text, it allows for a simpler and more direct 
analysis of log-enhanced contributions associated with this quantity. 

However, for the purpose of this appendix, we consider only observed regions which are not 
much larger than y (i.e. IR cutoffs not much smaller than the relevant momentum scale k). Under 
this assumption, Vc_ and Vc_ coincide and all that follows applies to both spectra. 

We can compare the two definitions of the power spectra as follows. Fourier expanding the two- 
point function of curvature perturbations, we can write 



d f k dk 



dlny 



2 ^U /"I 



d f sm{ky) p 



dlny J ky 



V^{k)dlnk . 



ko 



To pass from the first to the second line, we used the definition of the power spectrum in momentum 
space, given before eq. d?]). It is labeled by an F, to distinguish it from the analogous quantity in real 
space. In the third line, we made the IR cutoff explicit introducing a scale ko. We can expand V^ as 

/ k\ dV^ 
PfW=Pffe) + lu(-)ji^(y + .... (29) 

where kp is some pivot scale, the inverse of which is of the order of y: we can set y ~ fc~^. The second 
and higher terms are slow-roll suppressed with respect to the first term. Substituting this expansion 
in eq. (l28l) . we find 



koy 



^c(^) ~ ^c^(^p) — ', + slow- roll suppressed terms (30) 



10 



where the slow- roll suppressed terms are weighted by logarithms of order In {k^y). We choose k^ such 
that k^y ~ k^/kp <^ 1, but not as small as to lead to large logarithms in the slow-roll suppressed 
terms. This implies that the power spectra in coordinate and momentum space coincide, up to 
negligible slow-roll corrections: 

v^y) ^ v^a/y) ■ (31) 
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